arXiv: 1506.07119v3 [gr-qc] 15 Mar 2016 


A class of black holes in dRGT massive gravity and their 
thermodynamical properties 

Sushant G. Ghosh Lunchakorn Tannukij and Pitayuth Wongjun]^’® 

^Centre of Theoretical Physics, Jamia Millia Islamia, New Delhi 110025, India 
“^Astrophysics and Cosmology Research Unit, 

School of Mathematical Sciences, University of Kwazulu-Natal, 

Private Bag BfOOl, Durban fOOO, South Africa 
^Department of Physics, Faculty of Science, 

Mahidol University, Bangkok lOfOO, Thailand 
^The institute for fundamental study, 

Naresuan University, Phitsanulok 65000, Thailand 
® Thailand Center of Excellence in Physics, 

Ministry of Education, Bangkok lOfOO, Thailand 

Abstract 

We present an exact spherical black hole solution in de Rham, Gabadadze and Tolley (dRGT) 
massive gravity for a generic choice of the parameters in the theory, and also discuss the thermody¬ 
namical and phase structure of the black hole in both the grand canonical and canonical ensembles 
(for the charged case). It turns out that the dRGT black hole solution includes other known so¬ 
lutions to the Einstein field equations, such as the monopole-de Sitter-Schwarzschild solution with 
the coefficients of the third and fourth terms in the potential and the graviton mass in massive 
gravity naturally generates the cosmological constant and the global monopole term. Furthermore, 
we compute the mass, temperature and entropy of the dRGT black hole, and also perform thermo¬ 
dynamical stability analysis. It turns out that the presence of the graviton mass completely changes 
the black hole thermodynamics, and it can provide the Hawking-Page phase transition which also 
occurs for the charged black holes. Interestingly, the entropy of a black hole is barely affected and 
still obeys the standard area law. In particular, our results, in the limit mg —)• 0, reduced exactly 
to vis-a-vis the general relativity results. 
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I. INTRODUCTION 


The question of whether a mass term for the graviton field can be introduced existed 
in Einstein’s theory of general relativity since its inception. Massive gravity came into 
existence as a straightforward modification of general relativity by providing consistent in¬ 
teraction terms which are interpreted as a graviton mass. Such a theory can describe our 
Universe, which is currently undergoing accelerating expansion without introducing a bare 
cosmological constant. Massive gravity modifies gravity by weakening it at the large scale 
compared with general relativity, which allows the Universe to accelerate, while its predic¬ 
tions at small scales are the same as those in general relativity. Furthermore, if solution 
exists in this theory, it may also elucidate the dark energy problem. Hence, in recent years 
there were numerous developments in the massive gravity theories [IH6]. The first attempt 
was done, in 1939, by Fierz and Pauli jl]. They added the interaction terms at the linearized 
level of general relativity but later it was found that their theory suffered from discontinuity 
in predictions which was pointed out by van Dam, Veltman, and Zakharov, the so-called 
van Dam-Veltman-Zakharov (vDVZ) discontinuity [2H1]. This discontinuity problem in¬ 
voked further studies on the nonlinear generalization of Fierz-Pauli massive gravity. During 
the search for such generalizations, Vainshtein found the origin of the vDVZ discontinuity 
is that the prediction made by the linearized theory cannot be trusted inside some charac¬ 
teristic “Vainshtein” radius and he also proposed the mechanism for the nonlinear massive 
gravity which can be used to recover the predictions made by general relativity |5]. At the 
same time, Boulware and Deser found that such nonlinear generalizations usually generate 
an equation of motion which has a higher derivative term yielding a ghost instability in 
the theory, later called Boulware-Deser (BD) ghost [6|. However, these problems, arising 
in the construction of the massive gravity have been resolved in the last decade by first 
introducing Stiickelberg fields [7]. This permits a class of potential energies depending on 
the gravitational metric and an internal Minkowski metric. Furthermore, to avoid reap¬ 
pearance of the ghost in massive gravity, the set of allowed mass terms was confined and 
furnished perturbatively by de Rham, Gabadadze and Tolley (dRGT) [HllH]. They summed 
these terms and found three possibilities, viz. quadratic, cubic and quartic combinations of 
the mass terms. The dRGT massive gravity is constructed suitably so that the equations 
of motion have no higher derivative term, so that the ghost field is absent. However, these 
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nonlinear terms lead to complexity in the calculations and hence, in general, hnding exact 
solutions in this theory is strenuous. Nevertheless, recently, several interesting measures 
have been taken to obtain the spherically symmetric black holes in various massive gravities 
[Tnll26] . In particular, a spherically symmetric black hole solution with a Ricci flat horizon 
in four dimensional massive gravity with a negative cosmological constant was obtained by 
Vegh SOI. and was generalized to study the corresponding thermodynamical properties and 
phase transition structure IITHI3]. The spherically symmetric solutions for dRGT were also 
addressed in HSUS], the corresponding charged black hole solution was found in |Tn] and its 
bi-gravity extension, was found in na, which includes as particular cases as the previously 
known spherically symmetric black hole solutions. (See jlHH2Q], for reviews on black holes in 
massive gravity, see also IZU for the black hole solution in other classes of massive gravity). 

The main purpose of this paper is to present a new class of exact spherically symmetric 
black hole solutions including generalization to the charged case in dRGT massive gravity, 
and also to discuss their thermodynamical properties. It turns out that the solution discussed 
in this paper represents a generalization of the Schwarzschild solution that includes most 
of the known dRGT black hole solutions. The paper is structured as follows. In the next 
section, we review dRGT massive gravity. We present the modified equations of motion for 
dRGT massive gravity and a class of exact black hole solutions in section III. The calculation 
of the thermodynamical quantities associated with the dRGT black hole solution and the 
study of the phase structure of a black hole in the canonical ensemble approach are the 
main subject of section IV. The analyses in section IV are extended for the charged case in 
the section V, and Anally we summarize our results and evoke some perspectives to end the 
paper in section VI. We have used units which fix the speed of light and the gravitational 
constant via SttG = = 1. 

II. DRGT MASSIVE GRAVITY 

We begin by reviewing dRGT massive gravity, which is a well known nonlinear generaliza¬ 
tion of a massive gravity and is free of the BD ghost by incorporating higher order interaction 
terms into the Lagrangian. The dRGT Massive gravity can be represented as Einstein grav¬ 
ity interacting with a scalar held, and hence its action is the well-known Einstein-Hilbert 
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action plus suitable nonlinear interaction terms as given by j9j 


S 



R + mlU{g,r) , 


( 1 ) 


where R is the Ricci scalar and W is a potential for the graviton which modifies the gravita¬ 
tional sector with the parameter rUg interpreted as graviton mass. Moreover, the action is 
written in a unit such that the Newtonian gravitational constant is unity (thus, k'^ = Stt). 
The effective potential U in four-dimensional spacetime is given by 


W {a, <f) =U 2 + OL-M-i + ^4^4, 


( 2 ) 


in which and are dimensionless free parameters of the theory. The dependencies of 
the terms W 2 , W 3 and W 4 on the metric g and scalar fields 0 “ are defined as 


U, ^ []cr - [}C% (3) 

U3 = []Cf-3[IC][IC^]+2[IC% (4) 

W 4 = [JC]^ - 6[/C]2[/C2] + 8[}C][)C^] + 3[/C2]2 - 6[}C% (5) 

where 

= ( 6 ) 


where fab is a reference (or fiducial) metric and the rectangular brackets denote the traces, 
namely [/C] = and [JC^] = (/C”)(). The four scalar fields 0“ are the Stiickelberg scalars 
which are introduced to restore general covariance of the theory. Note that, generally, there 
are additional mass terms for the theory in higher-dimensional spacetime which are provided 
explicitly in One may recognize the interaction terms as symmetric polynomials of /C; 
for a particular order, each of the coefficients of possible combinations are chosen so that 
these terms will not excite higher derivative terms in the equations of motion. Actually, this 
definition of /C is not unique since it is possible to have the same action with a different 
definition of /C — the alternating action is given in 

To proceed further, we choose the unitary gauge 0“ = nni. In this gauge, the 

tensor gf^i, is the observable metric describing the hve degrees of freedom of the massive 
graviton. Note that since the Stiickelberg scalars transform according to the coordinate 
transformation, once the scalars are fixed, for example, due to choosing the unitary gauge, 
applying a coordinate transformation will break the gauge condition and then introduce 
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additional changes in the Stiickelberg scalars. Also, we redefine the two parameters 03 and 
Q !4 of the graviton potential in Eq. ([^ by introducing two new parameters ot and /?, as 
follows 


tts — 




1 — a 
12 


(7) 


By varying the action with respect to metric we obtain the modified Einstein field 
equations as 


— 0 , ( 8 ) 

where is the effective energy-momentum tensor obtained by varying the potential term 
with respect to 

+3,3 - ICICI + IKI" - l«'l} 

-^9,.{ra‘‘-3|C||K2| + 2|K»|}}. (9) 

In addition to the modified Einstein equations, one can obtain a constraint by using the 
Bianchi identities as follows 


= 0 , ( 10 ) 

where denotes the covariant derivative which is compatible with Henceforth, we 
shall use a and /?, instead of the parameters as and 04 . 


III. BLACK HOLE SOLUTION IN DRGT MASSIVE GRAVITY 


In this section, we will look for a static and spherically symmetric black hole solution of 
the modified Einstein equations ([^ with the physical metric ansatz 

( 11 ) 


ds^ = —n{r)dt^ + 2d(r)dtdr + —— + hir^dVl?^ 

fir) 


The solution is found and classified into two branches: d{r) = 0 or h{r) = h^r where ho is 
a constant in terms of the parameters a and (3 [271I29] . The most interesting branch is the 
diagonal branch, d{r) = 0, since it is simpler to analyze. For example, a class of a charged 
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black holes within the diagonal branch in dRGT massive gravity was also investigated in 

Ref. [30] . 

The exact solutions for this ansatz are complicated and thus it is difficult to use these 

solutions to analyze the properties of black hole. Note that one may simplify the solution 

by choosing some specihc relations of the parameters a and (3 for example a = —3/3 |16) . 
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Furthermore, a class of parameters satisfying (3 = ^ simply yields the Schwarzschild-de 
Sitter solution m (see also the analyses of such solution in Ref. j32f[35| i . 

It is important to note that most solutions are asymptotically de Sitter or anti-de Sitter. 
This is not surprising since at large scale the theory should recover the cosmological solution 
in which the graviton mass will play the role of cosmological constant to drive the late-time 
acceleration of the Universe. However, a class of the black hole solution in dRGT massive 
gravity (or in other classes of massive gravity, e.g. model in Ref. [36j) may encounter the 
issues of superluminality, the Gauchy problem, and strong coupling, (see j3Z] for the issues 
in dRGT and also 1381110] for those in another model of massive gravity) 

Since the hducial metric seems to play the role of a Lagrange multiplier to eliminate the 
BD ghost, one can choose an appropriate form to simplify the calculation. In the present 
work, we will follow [TOIITO] by choosing the hducial metric to be 

= diag(0, 0, c^, sin^ 9), (12) 


where c is a constant. With the choice of the hducial metric, the action remains hnite since 
it only contains non-negative powers of ffj,„ (see m. for more details). It is important to 
note that the ehective energy momentum tensor in Eq. (|^ is derived by requiring that the 


hducial metric must be non-degenerate. From Eq. (12), it is obvious that the hducial metric 
is degenerate and then one may not use the ehective energy momentum tensor expressed in 
Eq. (|^. However, as point out in |41| . one can use the Moore-Penrose pseudoinverse of the 
metric /C^ in order to hnd the ehective energy momentum tensor and it provides the same 
expression in Eq. 0, Therefore, one can use the ehective energy momentum tensor dehned 
in Eq. ([^ for the form of the hducial metric. Moreover, the results can be checked by using 
the mini-superspace action as usually done in cosmological analysis. One can obtain the 
equation of motion by using the Euler-Lagrangian equation for the variables, n and /. As a 
result, we found that the equations of motion still valid. 

For the physical metric, we will consider the diagonal branch of the physical metric by 
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setting d{r) = 0. In order to obtain a black hole solution, we will choose the function 
h{r) = r. Then, the physical metric can be written as 


ds = —n{r)dt + 


2 I I 2 7o2 

' + r ail . 


f(r) 


( 13 ) 


Up to this point, we have chosen to investigate just a class of the solutions to the dRGT 
massive gravity which possesses symmetries of our interest. Symmetries of solutions are of 
great important in massive gravity since they affect the number of degrees of freedom and 
the stability of the theory at times, which are also signihcant issues for many massive gravity 
models. As found in cosmological background, the number of degree of freedom crucially 
depends on the isotropy and homogeneity of the background physical metric as well as the 
form of the hducial metric Therefore, it is worthwhile to hnd, for this choice of the 

physical and the hducial metric,whether the number of degree of freedom still valids and 
each of them is healthy. This issue is out of scope of this work and we leave this investigation 
for further work. 


From the ansatz in Eq. (13), components of the Einstein tensor can be written as 

^ _ 1 


G r _ 


'f 7^2 ' 

f {rn' + n) 


nr‘2 


J.2 ’ 


' n' 

1 ^ 


( n” 

- - -|- 


— 

4n 

2r y 

l 2n 


n 


[n 


/32 


( 14 ) 

( 16 ) 

( 16 ) 


2nr 4n^ 

Computing the effective energy-momentum tensor in Eq. (§ with this ansatz, the tensor 


can be written as 


xl = - 


a{3r — c){r — c) 3l3(r — 3r — 2c 

n o 


x; = 


a(3r —c)(r —c) 3/?(r — c)^ 

o o 


r 

3r — 2c 
r 


XI = XI = 


q;(2c — 3r) ^ 3/?(c — r) ^ c — 3r 


( 17 ) 

( 18 ) 
( 19 ) 


Note that A* = A^. There are specihc values of the parameter c by which this effective 
energy-momentum tensor behaves like a cosmological constant. In particular, c = 0 simplihes 
each of the diagonal component of the effective energy-momentum tensor so that they depend 
only on the parameters of the theory, which are all equal constants. Actually, by setting 
c = 0 in Eq. ([^, the tensor /Cjj will be equal to the identity matrix leading to the fact 
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that the mass terms are all constants at the Lagrangian level which is corresponding to the 
cosmological constant term. This is a crucially different point between our model and one 
in Ref. m- In our model, the solution can be reduced to Schwarzschild-AdS/dS solution 
where the cosmological constant-like term can be expressed in terms of the graviton mass 
while the cosmological constant-like term in Ref. HH is introduced by hand and not related 
to the graviton mass. 

Substituting all components of Einstein tensor and effective energy momentum tensor 
into Eq. the modihed Einstein equations can be written explicitly as 


/' / 1 


a{3r — c){r — c) 3f3{r — 3r 

O O 


2 c 


/ (rn' -|- n) 


nr^ 


-^ 2=^9 


a{3r — c){r — c) 3(3{r — c) 


+ 


+ 


r 

3r — 2c 


n 


f I ^ 1 \ 

4n 2r / \ 2n 2nr 


[n 




= 


a(2c —3r) 3/3(c —r) c — 3r 


+ 


+ 


( 20 ) 

( 21 ) 


( 22 ) 


From Eq. (20), one can obtain the solution of / which can be written as 

2M A p 

/(r) = 1- + + 7r + C, 

r 3 


( 23 ) 


where 


A = 3m^ {1 + a + (3), 

(24a) 

7 = -cm^ {l + 2a + 3/3 ), 

(24b) 

C = c^m^ {a + 3/3 ), 

(24c) 


and M is an integration constant related to mass of the black hole. This solution incor¬ 
porates the cosmological constant term, namely A, naturally in terms of the graviton mass 
mg which should not be surprising since the graviton mass serves as the cosmological con¬ 
stant in the self-expanding cosmological solution in massive gravity. Moreover, this solution 
can be identihed with the known solutions in general relativity. Note that the solutions 
in alternative forms with other sets of parameters are shown explicitly in In the case 
mg = 0 we have the Schwarzschild solution, as expected. For c = 0 which sets 7 = ^ = 0, 
the solution can be classihed according to the values of a and (3. If (1 -|- a -1- /3) < 0, the 
















solution is in the form of Schwarzschild-de-Sitter while the case {1 + a + (3) > 0 yields the 


Schwarzschild-Anti-de Sitter solution. Using Eq. (20) and Eq. (21), one Ends that 


n'f = fn. 


(25) 


This equation implies that function / and n differ only by a constant. One can choose the 
constant to obtain a black hole solution such that 

n{r) = /(r) = 1 — + 'jr + (. (26) 

r 3 

It may be verihed by direct substitution this solution into Eq. ( [2^ . The dRGT solution out¬ 
lined here contains, for instance, the Schwarzschild solution [rrig = 0), the de Sitter/Anti-de 
Sitter solutions, the global monopole solution of general relativity, and thus also contains the 
monopole-de Sitter-Schwarzschild solution. Note that the last term, the constant potential 

corresponds to the global monopole term. A global monopole solution was introduced 
by Barriola and Vilenkin |53], and usually comes from a topological defect in high energy 
physics at early universe resulting from a gauge-symmetry breaking gniiiT]. However, in 
this solution, the global monopole is contributed from a graviton mass. In addition, we see 
that the solution here is similar to the four-dimensional solution found in Ref. m In Ref. 
m, the author obtained the black hole solution where the bare cosmological constant is 
presented while, on the other hand, our result incorporates the cosmological-constant-like 
behavior, namely the Ar^, in the black hole solution automatically due to introducing the 
potential term. 

However, this solution highly depends on the choice of the hducial metric; changing 
to other forms of the hducial metric will signihcantly affect the solution. This kind of 
dependency is one of the important properties of massive gravity. For example, from a 
cosmological point of view, one cannot have a nontrivial hat cosmological solution with a 
Minkowski hducial metric m, only the open FLRW solution is allowed 031 . where the 
FLRW solution with arbitrary geometry exists when the FLRW hducial metric is considered 
m- By generalizing the form of the hducial metric, the nontrivial cosmological solutions 
can be obtained |15] . 
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Figure 1. The graph in the left panel shows the domain region in (a,/3)-space for the black hole 
solution with three event horizons. The horizontal-blue shaded region corresponds to the region 
with nig = c = 1 and M = 2 while the vertical-black shaded region corresponds to the region with 
nig = c = 1 and M = 1. Note that there is an overlapping region where the black hole can have 
three event horizons for both values of the mass M. For the right panel, the graph shows the profile 
of /(r) (with M = 1). The red-solid curve corresponds to the case nig = 0, which recovers the 
Schwarzschild geometry while the blue-dashed curve corresponds to a set of parameters nig = c = 1, 
a = —3, and /3 = 2.1. The black-dotted curve corresponds to set of parameters as nig = c = 1, 
a = 10, j3 = 0.5. Together, these demonstrate the existence of three event horizons. 

IV. THERMODYNAMICS OF THE BLACK HOLE 


For black holes in de Sitter space, there exists more than one horizon, and the multi¬ 
ple horizons correspond to different thermodynamic systems. Next, we shall explore the 


thermodynamics of the dRGT massive gravity black hole solution given by the Eq. (23) 
by assuming that the black hole is a closed system, i.e., no particle (or charge) transfer or 
creation/annihilation. In other words, the black hole is assumed to be a canonical ensemble 
system. The horizons, if they exist, are given by zeros of = 0 |48| or 


2M A . 

1-h —r^ + jr + ( = 0, 

r 6 


(27) 


which may admit three real roots. For some classes of parameter setup, there may exist up 
to three horizons as shown Figure From the left panel in this figure, we use a region plot 


to hnd the region in (a,/3)-space for the existence of three positive real roots of Eq. (27) 
by setting nig = 1, c = 1. We also pick up two points in this region to show the prohle of 
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/(r) in the right panel of this hgure. We note that the gravitational mass of a black hole is 
determined by /(r_|_) = 0 , which in terms of the onter horizon radius r+ reads 

'■" A „ . . P8) 

(29) 


^ = y (^1 + + 7r+ + C ) , 

M f. 


— — (l + Tu^gCp' [(1 + Q; + /5)ry — (1 + 2 o^ + 3^)r+ + (o; + 3^)]) , 

C Zi 

where f+ = r+/c. For simplicity, one can work in dimensionless parameters with units of 
= 1. In this unit, the conditions for the positive value of the black hole mass can be 


written as 


a > 


(1 + /3)f2 - (1 + 3/3)f+ + (1 + 3/3) 


(r+ - 1)2 

/3 > — 1 and a is arbitrary for f+ = 1. 


for 7 ^ 1 , 


( 30 ) 

( 31 ) 


The Hawking temperature associated with the black hole is related with the surface gravity 
(k) via T = k/{2'k). The surface gravity in terms of the metric function reads k, = 
and hence the temperature T of the black hole becomes 

1 


T = 


47rr I 


(l + Ar^ + 27 r+ + C) 


( 32 ) 


Taking the limit mg = 0, one recovers the temperature for the Schwarzschild black hole 

T = 


dvrr 4 


( 33 ) 


The crucially different point in the temperature prohle compared with one obtained in 
the Schwarzschild solution is that it is possible to hnd the positive local minimum of the 
temperature as shown in the right panel of Figure In the left panel of this hgure, we 
use units of = 1 and show the region matching the requirement of positivity of the 
local minimum of the temperature, T+(min) > 0 , together with the positivity of the horizon 
size, r+ > 0. We also adopt a simple choice of the parameters in this region, such as 
(a = 1 , /3 = 0 . 2 ), to illustrate the existence of the positive local minimum of the temperature 
as shown in the right panel of FigureNote that, in the unit of = 1, the dimensionless 
version of the temperature and horizon size can be written as T = Tc and f+ = r_|_/c. These 
are the actual values of the quantities plotted in the right panel of Figure]^ For simplicity, 
we set c = 1 which leads to mg = 1 and T = T as well as f+ = r+. 

Next, we turn our attention to the important thermodynamic quantity associated with 
the black hole horizon which is its entropy (S'). The black hole is supposed to obey the hrst 
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Figure 2. In the left panel, the graph shows the domain regions in (a,/3)-space for the condi¬ 
tion > 0 corresponding to the vertically-red-striped region and for the condition r_|_ > 0 

corresponding to the horizontally-blue-striped region. The intersection corresponds to the region 
satisfying both conditions; > 0 and r_|_ > 0. We choose the point in the intersection region 


as a = l,/3 = 0.2 to illustrate the temperature profile of the black hole given by Eq. (32) in the 


T-^dM = j 


right panel. Note that we have used the unit of = 1 for both plots and also set M = 1 for the 
plot in the left panel. 

law of thermodynamics, or dM = TdS. To calculate the entropy, we use 

Integrating the above equation leads to 

S = vrr^. 

This can be written as S' = ^ with the area of the horizon = dvrr^, which is the 
famous area law. Interestingly, the graviton mass does not signihcantly affect the form of 
the entropy; it contributes only as a correction for the horizon radius which can be seen 


(34) 


(36) 


explicitly from Eq. (29). The behavior of the temperature and its corresponding horizon 


is shown in Figure in which there exists a local minimum of the temperature. This 
feature suggests that there should be a transition between two states; from the non-black 
hole or hot flat space state to a black hole. The transition was realized by Hawking and 
Page I5DI (see also ISl), who found the characteristics of the so-called Hawking-Page phase 
transition. The transition exists if it is thermodynamically spontaneous, or alternatively, 
globally thermodynamically stable, by evaluating the free energies between two states. In 
other words, this corresponds to evaluating the Euclidean actions of those two states where 
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the temperature is treated as a period of the imaginary time. Since in this case, there is no 
particle transfer, we compute the Helmholtz free energy 


F = M -TS, 

= h (l + 4’’+ + 7r+ + c) - X + ^^’'+ + d 

= + (36) 

The existence of globally thermodynamical stability of the black hole is determined by the 
condition F < 0. Therefore, the transition exists if 

Ar^ >3(1 + 0, (37) 

where the transition from the vacuum state to the black hole takes place at F = 0. 




0.0 0.5 1.0 1.5 2.0 

r+ 



Figure 3. This is the plot of the heat capacity of the black hole, corresponding to the black-solid 


line, given by Eq. (39) with the parameter setup nig = c = l,a = 1, (3 = 0.2. The heat capacity 
diverges at the radius to which the minimal temperature is associated. The graph also shows the 
temperature profile in dashed-black curve. 


In addition to the globally thermodynamical stability, one can determine the locally 
thermodynamical stability by examining the sign of the heat capacity. The heat capacity of 
the black hole is given by 


Fx = 


OM 

~w 


r=r+ 




r=r+ 


By using the relations in Eq. ( [2^ and Eq. (32) , the heat capacity becomes 

27rr^ (l + 27 r_,_ + Ar^ + Q 


C^ = 


Ar2 - (1 + C) 


Frl - (1 + C) 


(38) 


(39) 
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One can see from Figure that there is a particular horizon for which the correspond¬ 
ing heat capacity diverges. The divergence is due to the minimum temperature and the 
corresponding horizon is exactly the horizon of minimum temperature. Furthermore, re¬ 
quiring the locally thermodynamical stability of the black hole, the black hole must obey 
the condition 


Ar+> (1 + C). 


(40) 



Figure 4. The left panel is the region plot of (a, /3)-space which satisfies the stability conditions with 
given values of the parameters; nig = 1, c = 1, and M = 1. The blue-shaded horizontally-striped 
region corresponds to the values of a and (5 by which the condition for globally thermodynamical 


stability in Eq.(37) is satisfied while the vertically-red-striped region implies locally thermodynam¬ 


ical stability according to the Eq. (40). The intersection corresponds to the region which satisfies 
both globally and locally thermodynamical stabilities. For the right panel, the graph shows the re¬ 
gion satisfying both of the thermodynamical stability conditions in blue-shaded horizontally-striped 
region, the region satisfying the conditions > 0 and r+ > 0 in vertically-red-striped region, 

and the region in which three-event-horizon black hole exists in the skewed-black-striped region. 


depend on the parameters A and C, which are determined by the parameters of the massive 
gravity theory, namely, a and (3. We use the plot to hnd the allowed region in (a,/d)-space 
by using the unit of = 1 and setting M = 1. We show the validity of those parameters 
where both the stabilities are assumed in the left panel of Figure Note that one can 


express r+ in terms oi a, (3 and M from Eq. (29) and then substitute this expression into 
the stability conditions to hnd stability regions in {a,/3) space where M is held hxed to be 
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a positive constant. From this fignre, one can see that there exists an allowed region for the 
transition phase with the parameters a,f3 ^ 0(1). This snggests that massive gravity can 
natnrally provide the Hawking-Page phase transition withont reqniring fine-tnning of the 
parameters. For the right panel of Fignre]^ we combine three important regions according 
to the previons plots inclnding the regions satisfying thermodynamical stability conditions 
(intersection region in the left panel of Fignre |^, T+(inin) > 0 and r+ > 0 (intersection region 
in left panel of Fignre [^, and the existence of three horizons (region in the left panel of 
Fignre [^. From this fignre, one can see that the thermodynamically stable region of the 
black hole together with positive temperatnre and horizon size is not compatible with the 
region indicating the existence of three horizons. This is dne to the fact that the temperatnre 
is proportional to /'(r+) and the condition of existence of three horizons is that there exist 
two extremnm points snch that f'{r) = 0. This implies that there always exists a range of 
horizon (r_|_) at which the black hole temperatnre is negative if there exist three real black 
hole horizons. Therefore, we nse a set of parameters which give rise to one or two horizons 
to illnstrate the thermodynamical qnantities snch as temperatnre and heat capacity. 




a a 

Figure 5. In the left panel, the graph shows the stability region together with positivity of T+(min) > 
0 and r_|_ > 0 with different values of the parameter c. The vertically-red-striped region, the 
horizontally-blue-striped region, and the skewed-black-striped region correspond to the region with 
c = 0.8, c = 1 and c = 1.2 respectively. In the right panel, we choose a common value of parameter 
with j3 = 0.1 to illustrate the stability region in (q;,c) space. Both figures show that the greater 
the value of the parameter c, the smaller the region of stability. 

It is worthwhile to note that in the expression of dimensionless variables, snch as Eq. 
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(29), the dimensionless variable of r+ is f+ = r+/c. Therefore, the dimensionless horizon 
size is inversely proportional to the parameter c. This means that the greater the value of 
the horizon size, the smaller the value of the parameter c. In the left panel of Figure [^, 
we explicitly show the stability region in (a,/d)-space with different values of the parameter 
c such that c = 0.8, c = 1.0 and c = 1.2. Furthermore, we also show the allowed region 
in (a,c)-space by fixing (3 = 0.1 in the right panel. From this figure, it is found that the 
greater value of the parameter c corresponds to the smaller value of the horizon size and 
the smaller allowed region in the parameter space. Therefore, the phase transition tends to 
occur more easily at large horizon size. Note that in this case, the black hole is treated as a 
canonical ensemble system where particle transfer is prohibited. We will discuss the charged 
black hole case in the next section. 


V. CHARGED BLACK HOLE 


In this section, the black hole with non-zero charge and both the grand canonical aspect 
of the black hole, where the charge transfer is allowed, and the canonical ensemble point of 
view will be discussed. Hence, it will be interesting to consider the charged generalization 
of the above solution. The action in Eq. ([^ with Maxwell term reads 


S' = y 






(41) 


where is the Maxwell strength tensor and A^ is a vector poten¬ 

tial (the action is written in the Gaussian unit). We limit our study by considering a 
spherically-symmetric dRGT black hole filled only with a static charge which is accompa¬ 
nied by A^ = {A{r), 0, 0, 0). The presence of the static charge gives rise to a non-vanishing 
energy-momentum tensor as 

/A2 /A2 /A2 fA^' 


TH = diag 


Svrn’ Svrn’ Svrn’ Svm 


(42) 


Thus, one can rewrite the Einstein equations in Eq. (20) and Eq. (21) with the presence of 
a static charge as 


n 


/(r!t±d _ i = .rnlx: - 


nr^ 




n 


(43) 

(44) 
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where Xi and are given in Eq. (17) and Eq. (18). Taking these eqnations into acconnt, 


one can find that the constraint in Eq. (25) still holds for the presence of a non-zero charge. 
This constraint also simplifies the eqnation of motion of A{r). By reqniring the spherically 


symmetric solntion as in Eq. (13), one find the eqnation of motion of A{r) as 


1 ff 


n 


A" + A' - + - - - = 0 , 


2 V/ 


n 


A" + A- = 0, 
r 


(45) 

(46) 


where the constraint in Eq. (25) is applied in the calculation. This equation simply implies 
the solution of the form 


Air) = - + Eo, 
r 


(47) 


which is exactly an electrostatic potential in a generic electrodynamics, where k is an inte¬ 
gration constant. The corresponding electric field is 


^ k 

E{r) = -VA{r) = — 


(48) 


To determine the value of k, one may consider an electric field from a charge Q at large r 
where the spacetime is asymtotically fiat which, in the Gaussian unit, must take the form. 


E{r oo) ~ 


(49) 


Obviously, this implies that the integration constant k must be identified to the charge Q. 

), one finds a charged dRGT black hole solution as 
A 


By solving Eq. (43) and Eq. 

/(Q)(^) = 1 


2 M ^4 2 

-h ^ + —r -F yr -t- 


(50) 


where A, 7 , C, are defined similarly as in Eq. (24), note that we still have the relation in Eq. 
(26). In general, this solution may have up to 4 horizons, as illustrated in Figure]^ Again, 
the innermost 2 horizons are the horizons that can be found in the Reissner-Nordstrom 
black hole in the general relativity theory, where the others are the cosmological horizons 
associated with the existence of the graviton mass. 


Following the same procedures as in section pAd one can find the mass and the temperature 
evaluated at the black hole horizon as 

.- ( 51 ) 

(52) 


M = y ( 1 + — + -r+ + 7r+ + C 


T = 


47rr_i 


1 - % + Ar+ -h 27 r+ -F C 
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Figure 6. In the left panel, the graph shows the domain region in (a, /?)-space for the black hole 
having four event horizons. The horizontally-blue-striped region corresponds to the region with 
mg = 1, c = 1, Q = 0.2 and M = 2 while the vertically-black-striped region corresponds to the 
region with mg = 1, c = 1, Q = 0.2 and M = 1. For the right panel, the graph shows the 
profiles of /(r) with charge (the blue-dashed curve and black-dotted curve) compared with the 
Reissner-Nordstrdm solution (the red solid curve). The blue-dashed and black-dotted curves show 
an example of the existence of four horizons for the charged black hole, with parameter set mg = 1, 
Q = 0.2, c = 1, a = —3, /3 = 2.1, and mg = 1, Q = 0.2, c = 1, a = 10, /3 = 0.5 respectively. 

Since there is non-zero charge involved, the thermodynamics will be different. In detail, 
one can treat the black hole as an open system, i.e. a grand canonical ensemble, where the 
charge transfer is allowed while another can view the system as a closed one or a canonical 
ensemble where the charge is a non-zero constant. 


A. Grand Canonical Ensemble 


To treat it as a thermodynamical object, one can consider the black hole as a grand 
canonical ensemble system where the chemical potential is held fixed as p. = ^ (in the 
Gaussian unit). The corresponding temperature and entropy are 


T = _ 

^ 47rr_| 

S = Trri. 


(1 + ArJ. + 2jr+ + ( - , 


(53) 

(54) 


From Eq. (53), the temperature will be lower due to the effect of the chemical potential, p. 
In the grand canonical ensemble, the corresponding free energy is given as the Gibbs free 
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energy, 


G = M -TS - 

= lr+ (^l-'^rl + C- . (55) 

From this expression, one can see that the free energy shift to the smaller value due to the 
contribution from the chemical potential in the last term. In other words, the contribution 
from charge makes the free energy more negative. The sign of the free energy depends on the 
value of (l — + C ~ other words, there exists globally thermodynamical stability 

when 

Ar+ > 3 (1 + C - (56) 

is satisfied. Moreover, one can check locally thermodynamical stability by considering the 
heat capacity, 

_ 27rr2 [(l + 27 r+ + Ar^ + C) ^ SvrV^T 

Ar2 -(l + C) + /i2 - + ^ > 

Requiring locally thermodynamical stability, one should have 

Ar^ > (1 + C - . (58) 

For this grand canonical ensemble aspect, both of the stabilities depend not only on 
the parameters A and ( like the former analysis but also fi where the former two are again 
determined by a and (3. In the left panel of Figurej^, it shows the validity of those parameters 
where both the stabilities are assumed. From this figure, one may see that it is not difficult to 
find the allowed region with parameters a, /3 ~ 0(1) which suggests that the theory naturally 
provides the Hawking-Page phase transition. In the right panel of Figure a comparison 
of the allowed regions between charged and non-charged (/i = 0) cases is illustrated. From 
this plot, it is found that the existence of the charge makes the allowed region smaller. 
Furthermore, we can see the correspondence between the minimal temperature and the 
divergence of the heat capacity graphically in Figure 

B. Canonical Ensemble 

On the other hand, if charge transfer is prohibited, one can consider the black hole as 
a closed system, or a canonical ensemble with fixed non-zero charge Q. The mass and 
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a 


a 


Figure 7. The graph in the left panel is the plot of a and /3 which satisfy the stability conditions 
for the grand canonical ensemble with given values of the parameters; nig = 1, c=l, /r = 0.5, and 
M = 1. The horizontally-blue-striped region corresponds to the values of a and (3 by which both 


conditions for local stability in Eq. (58) and global stability according to the Eq. (56) are satisfied 
while the vertically-red-striped region satisfies the conditions Tg(min) > 0 r_|_ > 0. For the right 

panel, the vertically-red-striped region corresponds to the overlapping region from the left panel, 
in which both of the stability conditions and > 0, r+ > 0 are satisfied, compared with that 

obtained from the non-charged case, which is the blue-shaded horizontally-striped region. 




Figure 8. These are the plots of the temperature and the heat capacity of the charged black hole 
respectively where nig = 1, a = 1, /? = 0.2, c = 1, and n = 0.5. Both the temperature and the heat 
capacity profiles of the charged black hole are given by the blue thick curves compared with the 
black dashed curves corresponding to those of the neutral black hole. The divergence in the heat 
capacity plot corresponds to the minimal temperature of the black hole. 
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temperature are given readily by Eq. (51) and Eq. (52) respectively. Furthermore, the 
corresponding entropy still obeys the area law, 


S = vrr^. (59) 

Since the charge is fixed, the appropriate free energy in consideration is the Helmholtz free 
energy, which is 


F = M -TS, 

1 / A 2 , Q2\ 

= 4’'+ (1 - 3 rU C + 3 ^ j . ( 60 ) 

Similarly, the condition for globally thermodynamical stability is 

Arl >s(l + ( + 3 ^y ( 61 ) 

Furthermore, to examine the locally thermodynamical stability, the corresponding heat ca¬ 
pacity is computed. 


a = 


27rr2 (1 + Arl + 27^^+ + ^ _ SvrV^ T 


Ar^ - (1 + C) + 3 


Q2 


Arl - (1 + C) + 3 


<92 


(62) 


The condition for a locally stable black hole to exist is 

au > (1 + c - 3? 


(63) 

Here, the conditions for both stabilities depend on A and ( (or similarly, a and /3), as well 
as on the charge Q rather than the potential /i as in the grand canonical ensemble case, as 
shown in left panel of Figure Again, from the allowed region in this figure, the theory can 
provide the Hawking-Page phase transition naturally. From the right panel of this figure, 
it is also found that the existence of the charge makes the allowed region smaller, similar 


to the grand canonical case. This is implied by Eq. (52) where the region for which T > 0 
is reduced as the presence of charge. Similarly, the correspondence between the minimal 


temperature and the divergence of the heat capacity can be seen in Figure 10 


The effect of charge on the stability regions of both the grand canonical ensemble and 


the canonical ensemble compared with non-charged case are shown in Figure IT From this 
figure, one can see the that the effect of the chemical potential /i in the grand canonical 
ensemble and the charge Q in the canonical ensemble decreases size of the allowed region 


of the parameters. These can be seen from Eq. (52) and Eq. (53) since the contribution 
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Figure 9. The left panel is the region plot in (a, /?) space satisfying both of the stability conditions 
with given values of the parameters; rrig = 1, c = 1, Q = 0.2, and M = 1. The horizontally-blue- 
striped region corresponds to the values of a and (5 by which both conditions for local stability in 


Eq. (63) and global stability according to the Eq. (61) are satisfied while the vertically-red-striped 
region satisfies both of the conditions > 0 and r+ > 0. For the right panel, the skewed- 

black-striped illustrates the overlapping region in the left panel, satisfying both local and global 
stability conditions, Tc(inin) > 0 > 0 i compared with the one obtained in the non-charged 

case which is illustrated in the horizontally-blue-striped region. 


from the charge and chemical potential makes the positive temperature region smaller. In 
the canonical ensemble, we also found that the parameter region that satishes conditions 
for the existence of four horizons is not compatible with the stability region, similar to the 
non-charged case. The effect of the parameter c on the region of stability is also similar to 
the non-charged case. The stability region will increase where the horizon size increases or 
the parameter c decreases. 


VI. CONCLUDING REMARKS 

The dRGT massive gravity is a natural extension of Einstein’s theory of general relativity, 
providing mass to the graviton, and it is of a great arena for theoretical physics research. The 
dRGT massive gravity describes nonlinear interaction terms as a correction of the Einstein- 
Hilbert action and hence admits general relativity as a particular case. It is believed that 
dRGT massive gravity may provide a possible explanation for the accelerated expansion of 
the Universe that does not require any dark energy or cosmological constant. Hence, dRGT 
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Figure 10. These are the plots of the temperature and the heat capacity of the charged black hole 
respectively in the canonical ensemble view where rug = 1, a = 1, j5 = 0.2, c = 1, and Q = 0.2. 
Both the temperature and the heat capacity profiles of the charged black hole are given by the red 
thick curves compared with the black dashed curves corresponding to those of the neutral black 
hole. The divergence at large r_|_ in the heat capacity corresponds to the minimal temperature of 
the black hole. Moreover, the heat capacity diverges at two values of r_|_. 


massive gravity has received significant attention |52( including searches for black holes |18) . 
In this paper, we have obtained a class of black hole solutions in dRGT massive gravity, 
and studied the thermodynamics and phase structure of the black hole solutions. In the 
dRGT massive gravity outlined in this paper, there are three terms in the effective potential 
associated with the graviton mass. Furthermore, we note that due to the inclusion of the 
massive gravity term in the action, the Schwarzschild solution in the general relativity is 
modihed. Interestingly, it turns out that solutions to the Einstein held equations, such as 
the monopole-de Sitter-Schwarzschild become solutions in dRGT massive gravity for suitable 
choices of the parameters of the theory, where the coefficients for the third and fourth terms 
in the potential and the graviton mass in massive gravity naturally generate the cosmological 
constant and the global monopole term. The corresponding thermodynamical quantities are 
also changed. However, the black hole entropy is not affected signihcantly by the existence of 
the graviton mass, and still obeys the standard area law as in general relativity. Moreover, we 
consider the charged black hole solution in both the grand canonical and canonical ensembles 
to analyze the thermodynamics and phase transition. We have demonstrated through the 
calculation of the heat capacity and the free energy that there is a critical point where the 
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Figure 11. This plot shows the valid regions of the parameters a and /3 satisfying both global and 
local thermodynamical stability as well as the positive temperature condition and r_|_ > 0 where 
mg = c = 1, and M = 1. The vertically-red-striped region corresponds to the grand canonical 
aspect of the black hole with /r = 0.5, the horizontally-blue-striped region corresponds to the 
neutral black hole, and skewed-black-striped region corresponds to the canonical ensemble one with 

g = o. 2 . 


heat capacity diverges and a phase transition is possible without requiring hne-tuning of the 
parameters as shown in Figure]^ Even though it is possible to obtain three horizons in some 
region of parameter space, the region is still not compatible with the stability region. This 
implies that the phase transition will not occur when the black hole has three horizons. The 
presence of the charge will not change this argument, the phase transition will not occur 
when the black hole has four horizons. The presence of the charge affects the appearance 
of the Hawking-Page phase transition such that the allowed parameter region decreases as 


shown in Figure IT We also found that the phase transition tends to occur in the large 
horizon size in both the charged and non-charged cases. 


The black hole solutions obtained are immensely simplified due to the choice of the 
hducial metric as = diag(0, 0, c^, sin^ 0) and the choice of the Stiickelberg scalars. It 
will be interesting to apply the technique discussed here in other massive gravities to get 
black holes. It will also be interesting to consider the motion of particles in the background 
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of the dRGT massive black holes considered and to see how the graviton mass effects the 
eqnations of motion. These and related areas are for fntnre investigation. 
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Appendix A: Alternative form of dRGT action 

There is one of the alternative form of 4-dimensional dRGT massive gravity which is 
worthy of discnssion To consider the alternative form, it is nsefnl to write the ghost- 
free massive gravity action in a more general form. In an arbitrary nnmber of dimensions, 
n, the ghost-free interaction can be constrncted as 


n 



i=0 


where the term corresponds to the anti-symmetric contraction of the order of /Cjj as 


follow. 


Wo(/C) = 1, 
Wi(/C) ^/ccj. 




(A2a) 

(A2b) 

(A2c) 

(A2d) 

(A2e) 

(A2f) 

(A2g) 
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where the building block tensor is 


^u = \/9^''fabda(p'^du4)^. (A3) 

Note that here the square brackets acting on indices denote anti-symmetrization, 

[ , ] = ^{(even permutation of n indices) — (odd permutation of n indices)}. (A4) 

Generally speaking, this form of interaction can be expressed in an alternative form with a 
different definition of the building block tensor, namely 

n 

W(X) = Y^ciAm- (A5) 

j=0 

where 

X = \/ fabda(i)°-dy(t)^, (A6) 


For the 4-dimensional ghost-free massive gravity, the anti-symmetric contractions of the 
terms that are higher than 4*'^ order vanish which leaves the non-zero interaction term in 


Eq. (A2) as follows 


Wo(/C) = 1, (A7a) 

Wi(/C) = [/C], (A7b) 

W2(/C) = [/C]2-[/C2], (A7c) 

W3(/C) = [/Cj^^ - 3[/C][/C2] + 2[K% (A7d) 

U^{1C) = [Kf - 6[/C]2[/C2] + S[)C][}C^] + 3[/CY - G[}C% (A7e) 


To transform these terms to the different convention, introduced above, one can transform 
the coefficients Oj’s to Cj’s via the transformation matrix. 


(co) 


(l 

4 

12 

24 

24 

\ 


ao ^ 

Cl 


0 

-1 

-6 

-18 

-24 



tti 

C2 

= 

0 

0 

1 

6 

12 



tt2 

C3 


0 

0 

0 

2 

3 

2 

3 



«3 

VC4/ 


VO 

0 

0 

0 

1 

24 



V «4 / 


(AS) 


26 











or conversely via the inverse transformation matrix, 


do ^ 


/l 4 12 86 -1152 \ 


^ Co 



0 -1 -6 -27 720 


Cl 

02 

= 

0 0 19 -144 


C2 

03 


0 0 0 -| -24 


C3 

\ 04 / 


1^0 0 0 0 24 yi 


\C4/ 


(A9) 


In the case of 4-dimensional dRGT massive gravity, the zeroth order term; f/o, corresponds 
to the cosmological constant, then one can set Oq = 0 for simplicity. Moreover, the tadpole 
term Ui must vanishes to recover the Fierz-Pauli massive gravity in a linearized level [Tj, then 
= 0. Thus, the 4-dimensional dRGT massive gravity can be expressed in an alternative 
form as 


Cq — 12 (1 -|- 2q!3 -|- 204), 
C\ = —6 (1 -|- 803 -|- 4 tt 4 ), 
C2 = 1 -|- 60(3 -f 12cr4 , 

C3 = — 2 («3 + 014 ), 


C4 = 


0:4 


(AlOa) 

(AlOb) 

(AlOc) 

(AlOd) 

(AlOe) 


In addition, with the redefinition of parameters in Eq. Q, one can rewrite Eq. (AlO) as 


Co — 6 (1 -|- Cl -[-/?), 

(Alla) 

Cl = —2 (1 -|- 2ci -|- 3/3 ), 

(Allb) 

C 2 = 0 -|- 3/3 , 

(Allc) 

C 3 = ^ (1 - a - /?), 

6 

(Alld) 

^ 2 M • 

(Alle) 


Moreover, the black hole solution in Eq. ([23|) can be rewritten in this convention as 


where, alternatively. 


, 2M A 2 ^ 

/(r) = 1 - + -r^ + ir + 

r 8 


A = 


(Al 2 a) 


cm„ 


7 = 


-Cl, 


2 

C = c^mlc2. 


27 


(Al 2 b) 

(Al 2 c) 
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